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Abstract The Helioseismic and Magnetic Imager (HMI) and the Atmospheric 
Imaging Assembly (AIA) instruments onboard the Solar Dynamics Observatory 
satehite produce Doppler velocity and continuum intensity at 6173 A as weh as 
intensity maps at 1600 A and 1700 A, which can be used for hehoseismic studies 
at different heights in the solar photosphere. We perform a Hankel-Fourier anal- 
ysis in an annulus centered around sunspots or quiet-Sun regions, to estimate the 
change in power of waves crossing these regions of interest. We find that there 
is a dependence of power-reduction coefficients a on measurement height in the 
(-H I photosphere: Sunspots reduce the power of outgoing waves with frequencies v 

^ ^ lower than i> « 4.5 mHz at all heights, but enhance the power of acoustic waves 

in the range v w 4.5 — 5.5 mHz toward chromospheric heights, which is likely the 
signature of acoustic glories (halos). Maximum power reduction seems to occur 
near the continuum level and to decrease with altitude. Sunspots also impact 
the frequencies of outgoing waves in an altitude-dependent fashion. The quiet 
Sun is shown to behave like a strong power reducer for outgoing / and p-modcs 
at the continuum level, with a power reduction a ~ 15 — 20%, and like a weak 
power enhancer for p-modcs higher in the atmosphere. It is speculated that the 
surprising power reduction at the continuum level is related to granulation. In 
Dopplcr-vclocity data, and unlike in intensity data, the quiet Sun behaves like 
a strong power reducer for granular flows. 

Keywords: Instrument: SDO/HMI, SDO/AIA • Helioseismology 



1. Introduction 

The Helioseismic and Magnetic Imager instrument (HMI; Schou et ai, 2012) 
onboard the Solar Dynamics Observatory satellite (SDO) measures, in the Fe i 
line at 6173 A (e.g. Dravins, Lindegren, and Nordlund, 1981; Norton et ai, 2006), 
the motion of the solar photosphere. HMI samples the iron line at six wavelengths 
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symmetrical around the line center at rest. Line-of-sight (LOS) observables are 
produced every 45 s, which are Doppler velocity, LOS magnetic-field strength, 
Fe I line width, line depth, and continuum intensity. These are calculated from 
an algorithm based on the one used for the MDI instrument (Michelson Doppler 
Imager, Scherrer et al, 1995). A brief description of this algorithm is given in 
Couvidat et al. (2012). The Atmospheric Imaging Assembly instrument (AIA; 
Lemen et al., 2012) images the Sun in ten EUV and UV wavelength channels. 
It has conclusively been shown that two channels, at 1600 A and 1700 A, can 
be used for helioseismic studies {e.g. Hill et al., 2011; Howe et al., 2011, 2012; 
Rajaguru et al., 2012). 

The formation height of the signal in these two bands is « 430 and « 360 km, 
respectively, based on the study by Fossum and Carlsson (2005) for the TRACE 
mission. Images in these bands arc taken at a 24-s cadence. The formation height 
of the HMI Doppler velocity is « 150 km (Fleck, Couvidat, and Straus, 2011), 
while the formation height of the core of the Fe i line is « 302 km (Norton 
et al., 2006). Therefore, by combining HMI and AIA measurements, the entire 
photospheric layer and the lower chromosphere can be sampled (AIA UV bands 
are sensitive to heights up to « 800 km). This brings forth new opportunities to 
study this layer with the additional height information. 

In this paper, a Hankcl-Fourier transform (Braun, Duvall, and LaBontc. 1987) 
is applied to sunspot data, to measure the power reduction and/or enhancement 
in the / and p-oscillation modes crossing an active region. Braun, Duvall, and 
LaBonte (1987) established that sunspots "absorb" acoustic power, and this 
power reduction has been observed at different wavelengths, but no attempt has 
been made at measuring this reduction simultaneously at different atmospheric 
heights. Instead, different groups used individual datasets. For instance, Braun, 
Duvall. and LaBontc (1987, 1988) used Dopplcrgrams of the Fe i line at 8688 A 
formed in higher photospheric layers, Bogdan et al. (1993) used Dopplergrams 
of the Fe i line at 5576 A formed in the lower photosphere, Braun (1995) used 
Ca II K line intensity images formed in the chromosphere, and Zhang (1997) 
used intensity images at 3940 A sensitive to the middle and upper photosphere. 
Moreover, MDI and HMI continuum-intensity data appear not to have been used 
in any publication related to the Hankel-Fourier analysis. Even though Braun 
and Fan (1998) did use MDI intensities to measure meridional circulation, no 
power-reduction analysis was performed. 

The difference, if any, in power reduction at different heights should provide 
additional information regarding the physical mechanisms at play. Indeed, the 
complex interplay between magnetic field and plasma in the photosphere deeply 
impacts wave-propagation properties, probably in a height-dependent fashion. 
Whether the waves are partly running or completely standing, whether the 
mechanism responsible for power reduction operates at a given height or across 
a significant altitude or depth range, whether specific frequencies or wavelengths 
are affected at different heights, are but a few questions that can be investigated 
by Hankel-Fourier transforms using contemporaneous HMI and AIA data. 

Amongst all of the magnetic phenomena impacting acoustic-gravity waves, 
one of the earliest known is the reduction in amplitude of p-modes by active 
regions (Woods and Cram, 1981). Active regions also convert the acoustic and 
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surface gravity waves into magneto-acoustic-gravity (MAG) waves {e.g. Spruit 
and Bogdan, 1992; Gaily and Bogdan, 1993; Crouch et al., 2005; Khomenko 
and GoUados, 2006) that propagate higher in the atmosphere and deeper in 
subpliotospheric layers. Magnetic fields inhibit granulation, resulting in a lower 
acoustic-wave emissivity in sunspots {e.g. Parchevsky and Kosovichev, 2007; 
Hanasoge et ai, 2008), and lower the acoustic cut-off frequency {e.g. De Pontieu, 
Erdelyi, and James, 2004), allowing otherwise standing waves to travel upwards. 
Finally, magnetic fields modify acoustic-wave frequencies. 

Even though the photosphere is stably stratified, overshoot from the convec- 
tive zone produces a clear granulation and supergranulation signal in photo- 
spheric measurements that varies with altitude. Indeed, as pointed out by Howe 
et al. (2011), AIA 1600 and 1700 A intensity maps are much less contaminated 
by granulation noise than the HMI continuum. Therefore, noise properties are 
also height dependent, providing yet another rationale to performing helioseismic 
studies simultaneously on the different SDO observables. 

Several local-helioseismology techniques exist to study wave properties and 
infer the subphotospheric structure and dynamics {e.g. time-distance analysis, 
Duvall et ai, 1993; acoustic holography, Lindsey and Braun, 1998; ring-diagram 
analysis. Hill, 1988). These techniques complement or supersede the somewhat 
simpler Hankel-Fourier transform, recently fallen out of favor. Each one of these 
techniques provides a fragmentary picture of how strong magnetic fields in 
sunspots affect waves, and has been applied to data from various instruments, 
whose signals are not necessarily formed at the same height in the solar atmo- 
sphere. Here, the Hankel-Fourier transform is performed on signal from outside 
the sunspots {i.e. not contaminated by strong fields) to facilitate the interpre- 
tation of its results, and we use simultaneous measurements at four different 
altitudes. We study both sunspots and quiet Sun regions. In Section 1 we briefly 
remind the reader of the Hankel-Fourier analysis. In Section 2 we describe the 
data used. In Section 3 we present our results, and we conclude in Section 4. 

2. Hankel-Fourier Transform 

The Hankel-Fourier decomposition is extensively described in Braun, Duvall, 
and LaBonte (1987, 1988), and Braun (1995). Therefore, only a brief description 
will be given here. In polar coordinates (r, 9) and as a function of time t, the 
wave field (j) can be decomposed into components of the form 

^m{r,e,t) = e*('"''+"*)[A„(fc,c^)i/,W(fcr) +i3„(fc,^)i/,(?)(fcr)] (1) 

where m is the azimuthal order, k is the horizontal wavenumber, uj = 2t:i' is 
the cyclical frequency {v is the temporal frequency), and i/m'' and i/m'' are 
Hankel functions of the first and second kind, respectively. Goefficients Am and 
Bm are complex amplitudes of the radially ingoing and radially outgoing waves, 
respectively. As mentioned in Bogdan et al. (1993), Hankel functions are an 
approximation to the Legendre functions, valid for £ 3> m, where i is the angular 
degree (Braun, 1995). Here we only compute the power of ingoing and outgoing 
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waves, i.e. |^mp and |i?mP, and neglect the phases of these coefHcients. Indeed, 
phase measurement is a more comphcatcd matter and will be addressed at a 
later time. A measure, in an annulus centered on r = 0, of the reduction in 
power of outgoing waves am{k,uj), is obtained by 

am(fc,W =1--— n TT^ r 2 

where Nm{k, uj) is an estimate of the background power at k and uj for azimuthal 
order m. We are primarily interested in the power of oscillation modes. Therefore 
we subtract from the power spectra a background power resulting, mainly, from 
granulation and supergranulation. 

Quantity am{k,uj) is traditionally referred to as absorption coefficient in 
the literature. However, this coefficient includes the effects of not only power 
absorption, but also emissivity reduction arising inside sunspots (Chou et ai, 
2009b). The mechanisms responsible for these two effects are detailed in Chou 
et al. (2009a). Therefore, it is more appropriate to refer to am{k,u}) as a power- 
reduction coefficient (T.L. Duvall, Jr., 2012, private communication). 

To increase the signal-to-noise (S/N) ratio, an m-averaged outgoing- wave 
power-reduction coefficient a(/c,u;) is calculated, based on the sum of individual 
m power spectra (with or without background subtracted). For instance, with 
background subtracted, a{k,uj) can be defined as 

(Em \A„r{k,UjW) - N{k,Uj) 

where N{k,uj) is an estimate of the background calculated on the m-averaged 
power spectra. The same background is subtracted from ingoing and outgo- 
ing spectra. The calculation of a(fc,CL)) is somewhat sensitive to the way this 
background is removed, especially for the HMI continuum intensity where gran- 
ulation power is the highest. Therefore raw a{k,uj) values with no background 
subtraction are also presented here. Two techniques of background measurement 
are applied: one fitting a single background to the m-averaged power spec- 
tra (based on Braun, Duvall, and LaBonte, 1988), and one where a separate 
background power N„i is fitted to each individual m power spectrum (based on 
Braun, 1995). Both methods give similar results. In both cases, the logarithm 
of the power spectral density (PSD) is fitted by a 5th-order polynomial as a 
function of v, at each k. The / and pi to pg-mode ridges (oscillation modes with 
radial order n = to n = 9) are excluded from the fit. 

To further improve the S/N ratio, we may integrate the PSDs over w, or fc, 
or both, prior to calculating power-reduction coefficients. Thus, we define three 
quantities (here with the first method of background-power removal) , a{k), a{u}), 
and a, where 

^ 1 _ EJiErn\Brn{k,U;r)-N{k,U;)) 

^' EJiE„.\Amik,co)\^)~Nik,co))' 
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a{uj) ~ 1 



Y.k{{T.^\Am[k,w)V)-N{k,L0)y 



(5) 



and 



a 



1 



EfeEJ(E„JA„(fc,a.)P)-iV(fc,a.))- 



(6) 



A mask is usually applied to the PSDs, to only select power within the oscillation- 
mode ridges {e.g. Braun, Duvall, and LaBonte, 1988). Here, the mask only keeps 
the modes with n < 9 and above 2 mHz. 

The polar coordinate system is centered at a region of interest, e.g. the center 
of a sunspot, and the ingoing and out going- wave decomposition is calculated 
in an annulus with inner radius -Rmin selected to exclude this sunspot. For 
instance, Doppler-velocity measurements in strong magnetic fields are known 
to be unreliable [e.g. Rajaguru, Wachter, and Hasan, 2006), which makes it a 
better option to only use signal from outside strong-field areas. 

3. Data Used 

We use HMl continuum intensity, HMI LOS Doppler velocity, AIA 1700 A, and 
AIA 1600 A intensity maps. We also occasionally use HMI linc-corc intensity, 
formed by subtracting the Fc i line depth from the continuum intensity. The LOS 
magnetic flux in sunspots and quiet-Sun regions is estimated from HMI LOS 
magnetic- field strength (magnetogram). Each datacube is tracked for four days 
at the Carrington rotation rate using the mtrack software. A Postel projection 
is performed at each time step. The spatial resolution is dx = 0.09 heliocentric 
degrees, i.e. dx = 1.093 Mm (the HMI resolution is degraded, but this does not 
introduce significant aliasing in the power spectra), and the temporal cadence 
is dt = 45 s. A typical tracked cube has the following format: 384 x 384 x 7681, 
where the first two dimensions arc the spatial coordinates, and the last one is 
time. We select 62.5 h from each cube (most cubes have a few bad frames to be 
rejected). We set i?niin = 25.15 Mm, for all the sunspots to be contained inside 
the inner disk (except stmspot NOAA 11263, the largest studied which extends 
beyond -Rmin and is strongly asymmetrical). The outer radius of the annulus is 
Rmax — 209 Mm. There is a limit to the maximum azimuthal order m at which 
Am and Bm can be calculated: The inner-circle circumference 27ri?niin — 158 Mm 
corresponds to « 144 sampling points. Therefore, m (the number of nodes along 
the circle) should not exceed 144/2 = 72. Here we somewhat arbitrarily limit 
m to —46 < m < 46, because there is not much oscillation signal at higher m. 
The resolution in wavenumber (fc) is dk = 27r/(i?inax — -Rmin) = 0.0342 Mm~^, 
i.e. the resolution in angular degree I is equal to d£ — dk x Rq — 23.8 (where 
Rq ~ 696 Mm is the solar radius). The resolution in temporal frequency (v) is 
dv — 0.00444 mHz. Figure [T] shows the geometry used for the Hankel-Fourier 
analysis with respect to sunspot NOAA 11289. Table[T]lists all of the 15 sunspots 
studied. 

Figure [2] shows the PSDs (averaged over four quiet-Sun regions) for ingoing 
plus outgoing waves, as a function of £ and and for the HMI continuum, Doppler 
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Table 1. List of sunspots studied. 



NOAA number Start time (TAI) Stop time (TAI) 



11092 


2010, 


.07, 


.31. 


.12: 


:00: 


;00 


2010, 


.08, 


.04. 


.12: 


:00: 


:00 


11243 


2011, 


.07. 


.01. 


.20: 


:00: 


;00 


2011, 


.07, 


.05. 


.20: 


:00: 


:00 


11263 


2011, 


.08. 


.01. 


.12: 


:00: 


;00 


2011, 


.08, 


.05. 


.12: 


:00: 


:00 


11289 


2011, 


.09. 


.10. 


.00: 


:00: 


;00 


2011 


.09 


.14. 


.00: 


:00 


:00 


11306 


2011, 


.09. 


.29. 


.00: 


:00: 


;00 


2011, 


.10, 


.03. 


.00: 


:00: 


:00 


11314 


2011, 


.10. 


.13. 


.12: 


:00: 


:00 


2011, 


.10, 


.17. 


.12: 


:00: 


:00 


11317 


2011, 


.10. 


.14. 


.00: 


:00: 


:00 


2011, 


.10, 


.18. 


.00: 


:00: 


:00 


11352 


2011, 


.11. 


.19. 


.00: 


:00: 


:00 


2011, 


.11, 


.23. 


.00: 


:00: 


:00 


11384 


2011, 


.12. 


.24. 


.00: 


:00: 


:00 


2011, 


,12, 


.28. 


.00: 


:00: 


:00 


11388 


2011, 


.12. 


.31. 


.00: 


:00: 


:00 


2012, 


,01, 


.04. 


.00: 


:00: 


:00 


11408 


2012, 


.01. 


.25. 


.00: 


:00: 


:00 


2012, 


,01, 


.29. 


.00: 


:00: 


:00 


11410 


2012, 


.01. 


.30. 


.00: 


:00: 


:00 


2012, 


,02, 


.03. 


.00: 


:00: 


:00 


11419 


2012, 


.02. 


.16. 


.12: 


:00: 


:00 


2012, 


,02, 


.20. 


.12: 


:00: 


:00 


11420 


2012, 


.02. 


.16. 


.17: 


:00: 


:00 


2012, 


,02, 


.20. 


.17: 


:00: 


:00 


11423 


2012, 


.02. 


.27. 


.12: 


:00: 


:00 


2012, 


,03, 


.02. 


.12: 


:00: 


:00 



velocity, AIA 1700 A, and AIA 1600 A intensity data. The power spectra are 
clearly different. In particular, the high power in the granulation domain {v < 1.8 
mHz) on the HMI-continuum data is conspicuous, but it has already significantly 
dropped at the Dopplcr- velocity level. 

4. Results 

Four quiet-Sun regions are studied besides 15 sunspots. The quiet Sun was 
initially envisioned merely as a control of the Hankcl-Fourier transform code, 
but it turned out to present some new and interesting features on its own, and 
will be discussed after the sunspot results. 

4.1. Outgoing- Wave Power- Reduction and Enhancement in Sunspots 

To facilitate the comparison with past articles on power reduction in sunspots, 
we first calculate a{k) (Figure [21) and a{oj) (Figure H]). A background power is 
fitted to the m-averaged ingoing-plus-outgoing spectra, and removed. An / and 
p-mode mask selects the / and pi to pg ridges above 2 mHz for the computation 
of a{k) and a{u)). 

The different panels of Figure [3] look similar to the upper panel of Figure 
10 of Bogdan et al. (1993). Overall, sunspots do reduce the power of outgoing 
waves {a[k) > 0). There is an increase in a{k) with £, for £ < 500, and then a 
decrease for larger i. The peak power reduction occurs at the HMI continuum 
level and reaches 0.52. This is comparable to numbers quoted by Braun, Duvall, 
and LaBontc (1988) and larger than those of Bogdan et al. (1993). Any difference 
may be due to the sunspots themselves, to the way the background power is 
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Figure 1. Snapshot of a solar region centered on sunspot NOAA 11289 from the HMI con- 
tinuum intensity (upper left panel), absolute value of the HMI LOS magnetic-field strength 
(upper right panel), AIA 1700 A intensity (lower left panel), and AIA 1600 A intensity (lower 
right panel). The annulus inside which the Hankel- Fourier transform is calculated is defined 
by the two white circles. 



subtracted, to the / and p-niode mask applied, and to the higher m reached 
{e.g., Braun, Duvall, and LaBonte (1988) stopped at to = 5, Bogdan et al. 
(1993) stopped at m = 15, and Braun (1995) stopped at m = 20). 

The panels of Figure |3l characterizing four different heights in the photo- 
sphere, display different a{k) values, implying a dependence on height. The 
angular degree £ of peak power reduction slightly decreases with height from 
£ 500 at the continuum level to £ fa 400 higher in the photosphere (and even 
to lower £ for NOAA 11263 in AIA 1600 A data). 

On Figure 01 a{uj) peaks below i/ = 3 mHz, and decreases toward higher 
frequencies. The steepness of this decrease with ly, increases with atmospheric 
height for roughly v < 4.5 mHz. Moreover, a power reduction at the continuum 
level in the frequency range v w 4.5 — 5.5 mHz turns into a power enhance- 
ment {a(u)) < 0) in the upper photosphere. The existence of a negative a{uj) 
in this frequency range and for sunspot data has not been widely reported in 
previous publications. Braun (1995) mentioned a drop close to zero in power 
absorption near the photospheric acoustic cut-off frequency {ly w 5.3 mHz), 
and an increase at higher frequency, but no enhancement (although Figure 8 of 
Braun (1995) shows some amount of "emission" in the ranges £ = 329 — 390 and 
£ = All — 473, with only two sunspots analyzed). To the best of our knowledge, 
only Chen et al. (1996) explicitly mentioned a power enhancement, using K- 
linc data formed in the upper photosphere/lower chromosphere from the Taiwan 
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Figure 2. Maps of the power spectral density as a function of angular degree £ and temporal 
frequency u, averaged over four quiet-Sun regions in the HMI continuum intensity (upper left 
panel), HMI Dopplergram (upper right panel), AIA 1700 A intensity (lower left panel), and 
AIA 1600 A intensity (lower right panel). The color scale is truncated. 

Oscillation Network. With their annulus closest to the sunspot and similar to the 
one (i?min,fimax) used here, this power enhancement appears even at frequencies 
lower than in the present study, at least for one of their sunspots. 

The height dependence of a{u!) for waves with v < 5.3 mHz is difficult to 
explain if we assume that solar oscillations below the acoustic cut-off frequency 
(in the quiet Sun) are purely evanescent. However, propagation at frequencies 
below this cut-off was observed in sunspot penumbrae (Mcintosh and Jcfferies, 
2006; Rajaguru et al., 2007). Even in quiet Sun, Canfield and Musman (1973) 
early on discussed about low-frequency propagating waves. Roberts (1983) de- 
scribed how dissipation in the atmosphere lowers the cut-off frequency. Because 
radiative losses are important at the height of the HMI-continuum formation, 
there might not be a clear division between propagating and evanescent oscil- 
lations. More recently, De Pontieu, Erdelyi, and James (2004) investigated how 
inclined fields allow p-mode leakage, and Khomenko et al. (2008) showed how the 
5-min oscillations can propagate through the photosphere to the chromosphere 
in small vertical magnetic-flux tubes. Stodilka (2010) showed that granulation 
lowers the cut-off frequency in the solar photosphere, allowing for upward leakage 
of oscillation power. Finally, Howe et al. (2012) studied the phase of acoustic 
waves in and around active regions, using HMI and AIA observablcs. A phase 
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Figure 3. Power-reduction coefficient a{k) as a function of angular degree £ for 15 sunspots in 
the HMI continuum intensity (upper left panel), HMI Dopplergram (upper right panel), AIA 
1700 A intensity (lower left panel), and AIA 1600 A intensity (lower right panel). The orange 
curves are from a quiet-Sun region, the dashed lines are for sunspots with a LOS magnetic 
flux inside the disk of radius -Rmin larger than 8.4 X lO'^^ Mx, the dash-dotted lines are for a 
flux lower than 2.1 X 10^^ Mx, and the solid black lines are for fluxes in-between. Background 
power is removed and an / and p-mode mask is applied. The curves are smoothed. 

difference (implying propagation) in regions surrounding these sunspots is visible 
between HMI Dopplergram and AIA 1600 A intensity. All of these results, and 
the present paper, favor the idea that solar oscillations below the acoustic cut-off 
frequency are not purely evanescent. 

Figure [5] shows a as a function of the unsigned LOS magnetic flux inside 
the disk of radius i?inin, and for separate radial orders (n) of the oscillation 
modes (background power is removed), a is maximum at continuum level, and 
decreases with height: perhaps an indication that the mechanism responsible 
for power reduction operates mainly near, or below, continuum level. Maximum 
reduction occurs for pi and p2-modes, reduction is weaker for p3-modes, and 
minimum reduction occurs for /-modes. There is a dependence of a on magnetic 
flux: Larger fluxes produce greater reduction. This dependence seems stronger 
for p-modcs than for /-modes. Saturation might be reached for large fluxes, but 
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Figure 4. Power-reduction coefficient a(i/) as a function of temporal frequency i/ for 15 
sunspots in the HMI continuum intensity (upper left panel), HMI Dopplergram (upper right 
panel), AIA 1700 A intensity (lower left panel), and AIA 1600 A intensity (lower right panel). 
Background power is removed and an / and p-mode mask is applied. The curves are smoothed. 
The orange curves are for a quiet-Sun region. 



additional sunspots are needed to increase the sample size and access stronger 
fields. On FigurelHl the ratios of a in HMI Dopplergram to a in HMI continuum, 
and of a in AIA 1600 A to a in HMI continuum, both appear to increase with 
magnetic flux (although there is significant scatter in the data points). Thus, the 
region of outgoing- wave power reduction could extend higher in the photosphere 
as the flux increases. The main contributor to this power reduction is consid- 
ered to be mode conversion {e.g. Spruit and Bogdan, 1992; Gaily and Bogdan, 
1993; Crouch et al., 2005; Khomenko and Collados, 2006), where p-modes suffer 
absorption by partial conversion to MAG and Alfven waves. The slow MAG 
modes travel upward toward the chromosphere, along magnetic-field lines, while 
fast MAG modes travel upwards but arc progressively refracted and eventually 
reflected, probably in the lower chromosphere. Mode conversion occurs primarily 
in the equipartition layer, where Alfven speed equals sound speed. 
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Figure 5. Power-reduction coefficient a for the sunspots in the HMI continuum intensity 
(crosses), HMI Dopplergram (stars), AIA 1700 A intensity (diamonds), and AlA 1600 A 
intensity (triangles). Upper left panel is for the / modes, upper right panel is for the pi-modes, 
lower left panel is for the p2-rtiodes, and lower right panel is for the p3-modes. Background 
power is subtracted. 



Because some oscillation modes have their power enhanced by sunspots (Fig- 
ure [4]), the frequency ranges v < 4.5 mHz and 4.5 < v < 5.5 niHz are separated 
in Figure [T] This ensures that power enhancement for modes in the range 
v « 4.5 — 5.5 mHz does not impact power-reduction estimates at lower fre- 
quencies. The power enhancement is mainly present in AIA data (triangles and 
diamonds), i.e. the upper photosphere and lower chromosphere (although some 
sunspots also show a weaker enhancement in Doppler data), and seems to only 
have a weak dependence on the LOS magnetic flux. 

Due to the long duration of the tracked cubes, the high resolution and excel- 
lent quality of HMI and AIA data, it is possible to present the oscillation-power 
reduction/enhancement in sunspots in a more informative fashion by drawing 
2D maps of a{k,uj), as was done on Figure 5 of Bogdan et al. (1993) with 
lower resolution and a higher noise level. To increase the S/N ratio, the maps 
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Figure 6. Ratio of power reduction (a) in HMI Dopplcrgram to power reduction in continuum 
intensity (stars), and ratio of power reduction in AIA 1600 A intensity to power reduction in 
continuum intensity (triangle), as a function of the unsigned LOS magnetic flux of sunspots. 
Background power is subtracted and an / and p-mode mask is applied. Only the / and p-modes 
in the frequency range u < 4.5 mHz are selected. 
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Figure 7. Power-reduction coefficient (a) for the sunspots in the HMI continuum intensity 
(crosses), HMI Dopplcrgram (stars), AIA 1700 A intensity (diamonds), and AIA 1600 A 
intensity (triangles). Left panel is for the / and p-modes in the frequency range u < A.5 mHz, 
right panel is for the / and p-modes in the frequency range 4.5 < u < 5.5 mHz. Background 
power is subtracted and an / and p-mode mask is applied. 



of a{k,uj) for the 15 sunspots are averaged all together (despite each sunspot 
having a different size, shape, and magnetic flux). The result is shown in Figure 
IHl No background subtraction is performed. Therefore, a(fc, w) in the oscillation- 
mode ridges is underestimated. These maps present some striking features. First, 
the power enhancement mostly seen with AIA data on Figure S] and in the range 
v « 4.5 — 5.5 mHz is conspicuous. It is probably a manifestation of acoustic halos 
(Brown et at, 1992; Braun et ai, 1992). The term acoustic glory (Braun and 
Lindsey, 1999; Donea, Lindsey, and Braun, 2000) is also used when referring to 
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Figure 8. Maps of power-reduction coefficient a(k,uj) as a function of angular degree £ and 
temporal frequency u averaged over 15 sunspots in the HMI continuum intensity (upper left 
panel), HMI Dopplergram (upper right panel), AIA 1700 A intensity (lower left panel), and 
AIA 1600 A intensity (lower right panel). The color scale is truncated to —0.2 to +0.2. The 
maps arc smoothed to reduce the noise level. 

these extended halos where acoustic power surrounding a sunspot or a plage is 
enhanced. 

To try and determine where this enhancement occurs, Figure IHl shows the 
individual power-reduction coefficients (am) obtained for the AIA 1600 A data 
as a function of m, for / and p-modc ridges at all frequencies in the upper panel, 
and for p-modc ridges in the frequency range v = 4.5 — 5.5 mHz in the lower 
panel (with mask applied to select n < 9, but no background removed). Plotting 
function of azimuthal order m provides information about the spatial 
extent of the power-reduction/enhancement region. Indeed, as Braun, Duvall, 
and LaBonte (1988) point out, the "impact parameter" of a mode increases 
with to: The energy density of an incident plane wave falls off exponentially for 
a radial distance r < m/k. In other words, the higher m is, the less sensitive 
a mode is to the central part (small r) and the more sensitive it becomes to 
the peripheral part (large r) of the region under consideration. Therefore, in 
the upper panel of Figure [51 a drop in a„i with increasing to is consistent with 
the power-reduction region being contained inside the inner disk of radius i?min- 
am levels off at roughly m > 30 but does not reach zero for some sunspots, 
perhaps due to the presence of plages or other magnetized regions outside these 
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Figure 9. Power-reduction coefficient Om as a function of m for the sunspots in the AIA 
1600 A channel. Upper panel shows the power reduction integrated over all k and u (after 
applying the / and p-mode mask). Lower panel shows the power reduction in the frequency 
range v = 4.5 — 5.5 mHz (after applying the / and p-mode mask). The curves are smoothed 
to reduce the noise level. No background power is removed. 



sunspots: The region of power reduction seems to extend beyond the sunspots. 
Indeed, most sunspots are surrounded by plages visible in AIA data (for example, 
see Figure [ij, and plages are also present further inside the annulus where the 
Hankel- Fourier transform is performed (contaminating our results) . In the lower 
panel of Figure O (t^ = 4.5 — 5.5 mHz) the power enhancement occurs mostly 
at TO > 10. At lower m, the enhancement is less significant or even turns into 
a power reduction for some sunspots. This may imply that this enhancement 
occurs not toward the center, but rather toward the edges of — or even outside 
— sunspots. This supports the idea that it is a manifestation of acoustic glories. 
Moreover, this enhancement is not present at the continuum level (upper left 
panel of Figure [8|) , in agreement with the observation that there is no acoustic 
glories at this level (e.g. Hindman and Brown, 1998; Jain and Haber, 2002; Howe 
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et at, 2012). Acoustic glories are explained by Khomenko and CoUados (2009) as 
resulting from fast MAG waves being progressively refracted as they propagate 
upward in the sunspot. Because of this refraction, the waves escape the active 
region. Another possibility is that halos result from scattering of waves by the 
magnetic field of sunspots {e.g. Hanasoge, 2009). 

Figures [10] and [11] show cuts in the average PSDs of ingoing and outgoing 
waves (average over the 15 spectra from sunspot data), and as a function of ly 
(Figure [TUj) and £ (Figure [TT|) . In the lower panel of Figure [TU] (cut at £ ~ 357 
for the AIA 1600 A data), the frequencies (v) of outgoing p-modes &t v > 4.5 
mHz {jpQ to pg-modes) are reduced compared to ingoing modes. This frequency 
shift dv increases with v. A simple Gaussian fit of individual mode peaks returns 
a frequency difference (outgoing minus ingoing) of dv ~ —0.0098 mHz for pg, 
dv = —0.0209 mHz for py, and dv — —0.054 mHz for pg. 

Such an effect is absent (or less significant) at lower frequencies, and instead 
we obtain dv = +0.0095 mHz for pi-modes at ^ = 357. At ^ « 1400, dv of 
Pi-modes is close to zero. Finally, the /-modes seem to always show a small 
frequency increase {dv = 0.0104 mHz at £ = 357 and dv = 0.007 mHz at £ = 
1400) for outgoing waves. These results contradict Braun, Duvall, and LaBonte 
(1988) who observed no difference (within error bars) in the ridge positions of 
the ingoing and outgoing waves. The decrease in frequency of outgoing waves 
at higher v is absent at continuum or Doppler velocity level, and is visible only 
on AIA data. On the other hand, there could be a small positive dv at all 
frequencies in the HMI continuum (for pi at £ = 357, dv = 0.003 mHz, i.e. 
below the resolution in v). Frequency increase in AIA data for \ow-v modes 
is compatible with the results of Sun et al. (1997), who did observe such an 
increase on K-line images. They interpreted it as the effect of outgoing flows 
in and around sunspots. However, their Figure 1 seems to show a positive dv 
for V > 4.5 mHz, where we observe a negative one. The origin of the frequency 
decrease of higher-:/ modes in AIA data is more difficult to ascertain. 

This negative dv on AIA data partly contributes to the existence of a neg- 
ative a{uj) in the range v sa 4.5 — 5.5 mHz, but the lower panel of Figure 1101 
and the Gaussian fits of power peaks confirm that the power enhancement is 
genuine: The PSD is indeed larger for outgoing waves. Incidentally, the peak at 
V Ki ^.2 mHz highlights that a power enhancement of outgoing waves is present 
even in Doppler data. Therefore, this enhancement develops rather low in the 
photosphere and becomes stronger with height. 

Power reduction/enhancement could also partly be due to wave scattering by 
inhomogeneities in sunspots (Braun, Duvall, and LaBonte, 1987), and not only 
to actual power absorption/emission: Ingoing waves may be scattered into higher 
wavenumbers. Braun, Duvall, and LaBonte (1987) and Hanasoge (2009) consider 
such a change in wavelength at fixed frequency. Scattering of incoming waves by 
sunspots also affects their phases, e.g. Braun (1995). Figure [TT] shows cuts in the 
power spectra at v ^ b mHz and as a function of £. In AIA data, a Gaussian fit 
of individual mode peaks finds a very small angular-degree shift d£ of outgoing 
waves toward higher but the resolution in £ is not fine enough to confirm its 
existence {d£ = 1.934 at £ = 890 and ;/ = 5 mHz, and d£ = 2.23 at £ ^ 1179, 
well below the resolution in £). In any case, a small shift in wavelength and/or 
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Figure 10. Power spectral density of the ingoing (solid lines) and outgoing (dash-dotted lines) 
waves as a function of the temporal frequency v, averaged for the 15 sunspots, for continuum 
intensity data (upper panel), Dopplergrams (second row from top), AIA 1700A (third row 
from top), and AIA 1600 A (1 ower panel). These power spectra are cuts at £ = 357. The curves 
are smoothed. 
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Figure 11. Power spectral density of the ingoing (solid lines) and outgoing (dash-dotted 
lines) waves as a function of the angular degree i, averaged for the 15 sunspots, for continuum 
intensity data (upper panel), Dopplergrams (second row from top), AIA 1700A (third row 
from top), and AIA 1600 A (lower panel). These power spectra are cuts at u = 5 mHz. The 
curves arc smoothed. 
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frequency of outgoing waves is not the primary driver of the power enhancement 
observed in the upper photosphere in the v sa 4.5 — 5.5 mHz range, and of the 
power reduction observed at aU heights at lower frequencies. 

Noticeably, power- reduction coefficients calculated from Figure [TT] are higher 
(by w 0.01 at continuum level) than those shown on Figure [51 This difference 
arises because the average of the ratios outgoing PSD to ingoing PSD is not 
equal to the ratio of the average of outgoing PSDs to the average of ingoing 
PSDs. 

On Figure[T21 the 15 PSDs from sunspot data are averaged all together over k 
and normalized by their respective maximum values, to compare the frequency 
at which maximum power occurs. With AIA 1600 A the peak power of outgoing 
waves is shifted by dv = 0.087 mHz compared to ingoing waves (a Gaussian fit is 
used). At AIA 1700 A level, the shift is slightly smaller {diy = 0.082 mHz), and 
at Dopplcr-vclocity level it is only di' ~ 0.040 mHz. There is no significant shift 
at the continuum level. Therefore, dv is altitude dependent and results from the 
sunspots, since a similar plot with quiet-Sun data shows du ^ at all heights. 
Again, Sun et al. (1997) interpreted this overall frequency increase in the upper 
photosphere as resulting from outfiows in sunspots. It is also noteworthy that 
while the oscillation power peaks at v ^ 3.4 mHz with HMI continuum and 
velocity data, it peaks at v ^ 3.8 mHz with AIA data: probably reflecting the 
sensitivity to chromosphcric signal where the preferred oscillation period is 3 
min instead of 5 min at the photosphcric level. 

Figure [T2] shows function of m, for the /, pi, p2, and pa-modes 

separately, and at the HMI continuum level (only for the first four sunspots of 
Table [1] to make the figure more easily readable). a,„ decreases in all cases with 
m, but this decrease is steeper for p-modes than /-modes. Under the assumption 
that this dependence on m mainly reflects a difference in spatial location where 
power reduction occurs (but could also partly be due to, e.g., a lack of symmetry 
of the absorption region about the polar axis) , it appears that at the continuum 
level the outgoing-wave p-modc power is reduced mainly toward the sunspot 
center while the outgoing-wave /-mode power is reduced more evenly across the 
entire inner disk (and maybe even outside this disk). Indeed, /-modes show a 
smaller drop in am as m increases. For the p3-modcs, a™ levels off at m « 30. 

Figure [T3] is the same as Figure [T31 but at AIA 1600 A level. Here too, the 
decrease in a™ is steeper for p-modes than /-modes, and this power reduction 
turns into power enhancement for p-modes as m increases (only higher £ are 
accessible when m increases because the impact parameter of the waves increases: 
The power spectra become restricted to high-£ high-z/ modes). 

Finally, Figure [15] shows a as a function of the inner turning points Rq of 
modes, calculated in the ray-path approximation, for sunspot NOAA 11289. 
Only the first three radial orders n are plotted. Rq is defined as the depth at 
which c/i?o = where c = c(r) is the sound speed. Solar model S of J. 
Christensen-Dalsgaard provided c(r). There is a dependence of a on Rq that 
is similar for the three radial orders studied, a first increases with Rq, then it 
reaches a plateau, and it eventually decreases. This behaviour could be explained 
if the power reduction occurs in a layer extending underneath the photosphere. 
However, for n ~ S and at the continuum level, a increases until Rq « 14 Mm. At 



SOLA: ms.tex; 11 September 2012; 0:42; p. 18 



Oscillation Power in Sunspots and Quiet Sun from Hankel Analysis 




0.2t 

2.5 3.0 3.5 4.0 4.5 5.0 
V (mHz) 




2.5 3.0 3.5 4.0 4.5 5.0 

V (mHz) 




2.5 3.0 3.5 4.0 4.5 5.0 2.5 3.0 3.5 4.0 4.5 5.0 

V (mHz) V (mHz) 



Figure 12. Power spectral density of the ingoing (solid lines) and outgoing (dash-dotted 
lines) waves normalized by their maximum values, as a function of the temporal frequency v, 
averaged for the 15 sunspots, for continuum intensity data (upper left panel), Dopplergrams 
(upper right panel), AIA 1700 A (lower left panel), and AIA 1600 A (lower right panel). 



such a depth, the magnetic pressure is so low compared to gas pressure that it is 
unhkely that the field of sunspots has any impact. Therefore, the dependence of 
a on i?o may actually be a dependence on another parameter correlated with i?o : 
e.g., the angle of incidence of waves reaching their upper turning point varies 
with i?Oj and the inclination of p-modes with respect to magnetic- field lines 
impacts mode conversion. For AIA data, the power enhancement in each ridge 
occurs for the modes with the shallowest lower turning point, i.e. the modes 
closest to the surface. 

4.2. Outgoing- Wave Power-Reduction and Enhancement in Quiet Sun 

On Figures [3] and HJ a{k) and a(a;) are plotted for a quiet-Sun region (in orange). 
HMI Dopplergrams and AIA data show a weak power enhancement peaking 
around I ~ 600, as observed first by Bogdan et al. (1993) and later confirmed 
by Braun (1995) and Chen et al. (1996), and a weak power reduction at higher 
£. Bogdan et al. (1993) explained this enhancement by assuming that there is a 
statistical equilibrium at the solar surface where power absorption and emission 
are balanced. Modes that are truly global in nature have their power reduced 
when crossing a sunspot or a pore somewhere on the solar surface, and this 
will create a power emission at some other location in the quiet Sun where the 
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Figure 13. Power-reduction coefficient am for four sunspots in the HMI continuum intensity 
as a function of the azimuthal order m. Background power is not subtracted. The upper left 
panel is for the /-modes only, the upper right panel is for the pi-modes, the lower left panel 
is for the p2-iiiodes, and the lower right panel is for the p3-modes. The dashed lines are a 
second-order polynomial fit to the data. 

Hankel-Fourier transform is performed. Chen et al. (1996) also mentioned the 
existence of a center-to-hmb effect, an excess of power for waves propagating 
from Sun center toward the hmb, contributing to a negative a. 

However, HMI continuum data on Figure [3] show a different picture. There 
is a nearly constant and positive a{k) in the range ^ k, 200 — 1500, with an 
amplitude of a{k) w 15% (when background power is subtracted). This some- 
what puzzling result, opposite to what happens higher in the photosphere, was 
confirmed by running the Hankel-Fourier code on three other quiet-Sun regions 
observed at different dates and latitudes. We also ran a different code provided 
by Douglas C. Braun (2012, private communication), and performed the Hankel- 
Fourier analysis on a tracked cube of full-disk MDI continuum-intensity data to 
confirm that HMI and MDI continuums behave consistently. Finally, as a test, 
we processed HMI intensity images prior to the Hankel transform analysis: The 
limb darkening was (mostly) removed and a running difference of the images was 
taken (to get rid of low- frequency granulation signal). This pre-processing had 
only a minor impact on a. 

Figure [16] shows the outgoing-wave power reduction a(fc,a;) maps averaged 
over four quiet-Sun regions. The continuum-intensity map stands out as being 
the only one showing a systematic power reduction in the / and p-mode ridges. 
Other maps show a residual reduction in the power of /-modes, but no reduction 
and even a weak power enhancement in some p-mode ridges. The source of 



SOLA: ms.tex; 11 September 2012; 0:42; p. 20 



Oscillation Power in Sunspots and Quiet Sun from Hankel Analysis 




Figure 14. Power-reduction coefficient am for four sunspots in the AIA 1600 A intensity as a 
function of the azimuthal order m. Background power is not subtracted. The upper left panel 
is for the /-modes only, the upper right panel is for the pi-modes, the lower left panel is for the 
P2-modcs, and the lower right panel is for the p3-modes. The dashed lines are a second-order 
polynomial fit to the data. 

strong power reduction at the continuum level is unknown, but it is suggested 
by Thomas L. Duvall (2012, private communication) that granulation plays a 
role. Indeed, Hankcl-Fourier analysis performed on HMI line-core intensities of 
the four quiet-Sun regions shows no significant power reduction in the p-mode 
range. The main difference between continuum and line-core intensities is a « 200 
km shift in the height of formation, resulting in significantly less granulation 
signal at the line-core level. The existence of this ubiquitous power reduction at 
the continuum level seems to invalidate the assumption in Bogdan et at (1993) 
regarding a statistical balance between emission and absorption at the solar 
surface, when height is not taken into account. 

Kostyk and Shchukina (1999), Khomcnko et al. (2001), and Kostyk, Shchuk- 
ina, and Khomenko (2006) showed that the amplitude of the 5-min oscillations 
is much larger over intergranulcs than over granules and that oscillations are 
suppressed over granules. Roth et al. (2010) confirmed that waves excited in the 
intergranular lanes lose energy while probing a granule. Kostyk, Shchukina, and 
Khomenko (2006) found that even though there are waves propagating upwards 
over both granules and supergranules, the waves over granules are closer to 
standing waves than those over intergranular lanes. These results might partly 
explain the positive a(/c,cj) obtained with the Hankel-Fourier analysis. Indeed, 
numerous granules are present inside the inner disk studied. With a typical 
diameter of 1 Mm, w 630 granules are expected at any time in the disk of radius 
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Figure 15. Power-reduction coefficient a for sunspot NOAA 11289 as a function of inner 
turning point of the pi (leftmost curves), p2, and p3 (rightmost curves) modes. Background 
power is not subtracted. Shown are the HMI continuum intensity (upper left panel), HMI 
Dopplergram (upper right panel), AIA 1700 A (lower left panel), and AIA 1600 A (1 ower right 
panel). The curves are smoothed. 



i?min- Acoustic wavcs excited outside this disk and traveling toward its center 
might have their power reduced when crossing the central part of those granules. 
However, this does not explain why this power reduction is asymmetrical between 
ingoing and outgoing waves. 

A plot of a as a function of mode lower turning point (Rq) shows no de- 
pendence (or only a weak one) on Rq for Dopplcr velocity and AIA data, as 
expected. However, at the continuum level, a increases with Rq up to Rq ?» 6.5 
Mm for n = 3 (the data points are very scattered). This might indicate that the 
mechanism responsible for power reduction at the continuum level operates in a 
shallow layer underneath the solar surface, maybe related to the typical height 
of a granule. 
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Figure 16. Maps of power-reduction coefficient a(k,uj) as a function of angular degree £ and 
temporal frequency u, averaged over four quiet-Sun regions in the HMI continuum intensity 
(upper left panel), HMI Dopplergram (upper right panel), AIA 1700 A (lower left panel), and 
AIA 1600 A (1 ower right panel). The color scale is truncated to —0.2 to +0.2. 



Another conspicuous feature of Figure [12] is the strong power reduction below 
the /-mode ridge on HMI Dopplergrams, absent, or weak, on intensity maps. 
This power reduction reaches maximum below v = 2 mHz and peaks at the 
angular degree i w 950 (wavelength of w 4.5 Mm) corresponding to granular 
flows {e.g. Hathaway et ai, 2000). Therefore, in Doppler signal and unlike in in- 
tensity, the quiet Sun reduces the power of outgoing granular flows. The reason is 
currently unknown, but is probably somehow related to wave properties specific 
to Doppler measurements, vs. intensity measurements. This power reduction 
is stronger in quiet Sun than in sunspots (see Figure [5]) , where granulation is 
inhibited by strong magnetic fields. 

Finally, it was suggested that we run the Hankel-Fourier analysis on nu- 
merical simulations of the upper convective zone (T.L. Duvall, Jr., 2012, pri- 
vate communication) to try and understand the origin of the power reduction 
in / and p-modes in the continuum data, and in the granular-fiow domain 
in Doppler data. To this end, the 96 Mm-wide MHD simulations of Robert 
Stein were used {e.g., Stein and Nordlund, 2000; Benson, Stein, and Nordlund, 
2006): Data cubes of intensity at the solar surface can be downloaded from 
,http:/ /st einr.pa.msu.edu/^bob/data.html, with a cadence of dt = 60 s and a 
spatial resolution of dx « 0.095 Mm. About 8 h of data are available. Partly 
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because of the lower duration of the dataset, and mostly because of its relatively 
small size (96 x 96 Mm^), the Hankcl- Fourier coefficients are noisy. There does 
not seem to be any power reduction for outgoing oscillation modes present in the 
simulation, at least not above noise level (rms variation for a{k,Lo) is w 0.18), 
despite the conspicuous presence of granulation signal. Similarly, the Hankel- 
Fourier analysis run on simulated Doppler- velocity data finds no outgoing-wave 
power reduction in the granulation domain (or it is below noise level: The rms 
variation for a(fc,cj) in Doppler data is also ~ 0.18). 

5. Conclusion 

High temporal-cadence and high-resolution data taken simultaneously at differ- 
ent wavelengths by the HMI and AIA instruments onboard the SDO satellite 
make it possible to seismically probe the lower layers of the solar atmosphere 
and to study wave propagation and properties from the continuum level to the 
lower chromosphere. A Hankel-Fourier analysis was applied to 15 sunspots and 
four quiet-Sun regions observed by SDO. As was already known, sunspots reduce 
the power of oscillation modes crossing them at all heights in the atmosphere for 
frequencies lower than i/ « 4.5 mHz. This is likely a result of mode conversion by 
magnetic fields, reduced emissivity in the sunspots, mode scattering, etc. This 
power reduction increases with the LOS magnetic flux. The outgoing-wave power 
reduction seems to be stronger at the continuum level than higher in the atmo- 
sphere (when background power is subtracted prior to computing a), perhaps 
indicating that the source of reduction is located low in — or even below — 
the photosphere. The height dependence of the out going- wave power-reduction 
coefficient (a) also favors the idea that part of the solar oscillations below the 
acoustic cut-off frequency (in the quiet Sun) do propagate upward. On sunspot 
data, we observe a power enhancement higher in the photosphere for waves in 
the frequency range v « 4.5 — 5.5 mHz. This enhancement is barely visible in 
Doppler-velocity data, but is strong at the AIA 1700 and 1600 A level. Such a 
negative a was reported by Chen et al. (1996). This is likely the signature of 
acoustic glories observed around sunspots. A weaker power enhancement is also 
present in quiet-Sun data of the upper photosphere, as was reported by Bogdan et 
al. (1993) and Chen et al. (1996). Sunspots also shift the frequencies of outgoing 
waves, depending on the altitude and the angular degrees and frequencies of the 
modes. The frequency increase observed for some modes in AIA data is similar to 
the results of Sun et al. (1997), interpreted as evidence of outflows in sunspots. 

In the quiet Sun, there is a striking dependence of the power-reduction co- 
efficient on height. At the continuum level, the quiet Sun behaves like a power 
reducer for p and /-oscillation modes. This puzzling behaviour, never reported 
before as far as we know, is suggested by Thomas L. Duvall, Jr. (2012, private 
communication), to be somehow related to the presence of granulation. Indeed, 
granulation power is prominent in HMI continuum data but quickly vanishes 
higher in the atmosphere. Closer to the chromosphere, the quiet Sun becomes a 
power enhancer in most p-mode bands, while continuing to absorb /-modes. In 
Doppler-velocity signal, there is a strong outgoing-wave power reduction in the 
domain of granular flows. A weaker reduction is also present in sunspot data. 
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Hankel-Fourier analysis performed on Doppler velocity and continuum data 
from one of R. Stein's upper-convective-zone MHD simulation shows no power 
reduction (above noise level) at the continuum level in the / and p-mode bands, 
and no power reduction in Doppler data at the surface in the granular-flow 
domain. These simulations do have a clear granulation signal though, proving 
that the relations between granulation and the observed power reduction in the 
HMI continuum, and granulation and the power reduction in granular flows at 
Doppler- velocity level, are not straightforward. 

This mostly observational study leaves many questions unanswered. Future 
work will focus on the origin of the power reduction of / and p-modes in the 
quiet Sun at the continuum level (and the power enhancement higher), and on 
the reduction in granulation power in Doppler data. We also plan to compute 
the phase differences between ingoing and outgoing waves, to gather information 
about their running or standing nature, and about scattering phenomena. More 
sunspots need to be analyzed, to access larger magnetic fluxes. Hankel-Fourier 
analysis should also be performed on line-core intensity data in a systematic way, 
to access yet another height in the photosphere, and to emphasize the difference 
between intensity and Doppler-velocity signals (as line-core and Doppler sig- 
nals are formed at similar altitudes). Several seismic studies have recently been 
performed combining HMI and AIA data, and a comparison of their results is 
warranted: this could help interpret some of the results obtained here. 
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